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ABSTRACT 

Context. Calculation of stellar models with varying degrees of d-element enhancement. 

Aims. We investigate the influence of both a new generation of low-temperature opacities and of various amounts of a- 
element enhancements on stellar evolution models. 

Methods. New stellar models with two different a-element mixtures and two sets of appropriate opacity tables are com- 
puted and compared. The influence of the different mixtures as well as that of the improved generation of opacity tables 
is investigated. 

Results. It is found that around solar metallicity the new opacity tables have a drastic influence on stellar temperatures, 
which is mainly an effect of the new low-temperature tables, and not of variations in a-element enhancement factors. The 
latter, however, influence stellar lifetimes via systematic opacity effects at core temperatures. We trace the reason for the 
low-temperature table changes to errors in the old tables. 

Conclusions. Variations in a-element abundance ratios affect the main-sequence properties of super-solar metallicity 
stars significantly. Red giant branch effective temperatures depend only slightly on the specific mixture. Our older low- 
temperature opacity tables were shown to be erroneous and should no longer be used for stellar models with near- or 
super-solar metallicity. Corrected tables have already been produced. 

Key words. Stars: interiors - evolution - abundances - low-mass 
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1. Introduction 

It is now widely accepted that the chemical composition of stars can be classified by two main classes of 
metal mixtures: one with nearly solar metal ratios and another one, in which the so-called a-elements are 
enriched, while most other elements, in particular C, N, and those of the iron-group, are present in near-solar 
element ratios. The a-elements are those produced by or-captures on oxygen in Supernova II explosions; they 
are O, Ne, Mg, Ca, Si, Ti. These elements are enhanced over their solar relative number fractions by factors 
of 1.5-3 (or 0.2-0.5 dex on a logarithmic scale). Since oxygen is contributing around 50% of the metals, 
its a-enhancement can significantly increase the total metallicity. Using one element, in particular iron, as 
the metal indicator, can therefore be misleading in evaluating the total metallicity of stars. In the traditional 
picture of galactic populations. Pop I is characterized by a solar metal mixture, while Pop II shows significant 
a-enhancement. However, the correlation with total metallicity or [Fe/H] is not at all a strict one, since the 
relative a-element abundances depend on the star formation history and rate. It is known that in elliptical 
galaxies and possibly the galactic bulge stellar compositions with solar or super-solar iron abundances and 
strong a-enhancement exist, indicating intensive and short-lived star formation histories. F or a discussion 
of g-e lement evolution and its connection with galacto-chemical evolution, see the review bv lWheeler et alJ 
(Il989l) . 

Calculation of stellar models therefore requires for each value of total metallicity (Z) at least these two 
basic internal metallicity mixtures: a solar-scaled and an a-enhanced one. This requirement has been realized 
in the early nineties and since then the number of stellar libraries which also take into account a-enhancement 
has steadily increased. While the inclusion of individual element abundances is trivial for the nuclear energy 
generation, it is non-trivial for an elaborate equation of state and for opacities. In the former case, one trusts 
that as long as hydrogen and helium are dominating the gas, the individual element abundances within the 
metals is of negligible influence. However, for the opacities this no longer holds true since the various ele- 
ments contribute very differently to the total absorption. As the calculation of opacities has developed into a 
specialised branch, this requires separate calculations by groups active in this field or flexible programs for 
end-users. 

In the past we have tried to include a-enhancement in the opacity tables we used as best as possi- 
ble. In iBencivenni et al an a-enriched mixture was first taken into account using the Los Alamos 
Astrophysical opa city library ( Huebne r et al. 1977) for temperatures above 1 eV. Similar tables were used in 
IWeissetalJ ll l995h . where for the first time a-enhanced stellar models for elliptical galaxies were calculated 
with almost consistent input physics. However, at low temperature the opacities were still solar-scaled. While 
it was shown by Salaris et al. ( 1993) that at low metallicities opacities with solar-scaled abundances and es- 
sentially the same total metallicity are presumably a good replacement for a-enhanced tables, this could not be 
assumed for near-solar or even super-solar metallicities. The final step became possible when low-temperature 
opac ity tables, which include molecular absorption processe,s were made available by Alexander & Ferg uson! 
(11994 AF94), both for solar and a-enhanced metal mixt ures. These could then b e combined with high- 
temperature tables for completely identical mixtures from ('iglesias & Ro2ers"l996', OPAL; IR96). The en- 
hancement for the various a-elements are variable in this metal mixture (Salaris & Weiss 1998). These tables 
were used by the authors for globular cluster calculations dSalaris & WeisJl998l) demonstrating also that the 
approximation bv lSalaris et all (1 19931 ) begins to fail around [Fe/H] = -0.7, the iron content of 47 Tuc. A 
complete set of models with metallicities rangingjrom Z = 0.008 to 0.070 using the same opacity tables was 
added to the Padova stellar library bv lSalasnich et alJ (l200(j) . Similar, but independent o pacity tables from 
the same sources were used, e.g., in the models by the Victoria and Yale-Yongsei groups dVandenBer^ et alJ 
boOO: Yi et al. 2003). Here, all a-elements have been enhanced by the same factor, ranging in the different 

sets from 0.3 to 0.6 dex. 

As demonstrated first by I Salaris et al.l (Il993l) . and verified by, e.g., ^ Salaris & WeissI (Il998h and 

IVandenBerg et al. (^00), to first order a-enhanced opacities at low metallicity and for low mass stel- 
lar models can be replaced by scaled-solar ones for the same metal fraction Z, provided that the ra- 
tio (Xc -H Xn -H Xo H- XNe)/(XMg + Xsi -H Xs H" Xca + Xpe) - whcrc Xi dcnotcs the mass fraction of el- 
ement i - is approximately the same in both mixtures (in spectroscopic notation this corresponds to 
[{Xc + Xf.) + Xo + XNe)/(XMg + Xsi + Xs + Xca + Xps)]~ 0). As a consequence, differences in the internal 
distribution of a-elements are of minor significance, as long as this condition is satisfied. The extension of 
this result to solar or super-solar metallicities has however never been investigated. 

For a population synthesis project (CoeUio et al. 2006, in preparation), we decided to compare stellar 
models of low-mass stars using the previously introduced a-enhanced tables with varying enhancement factors 
to other models with new opacity tables, in which all a-elements are enhanced by a constant factor of 0.4 dex. 
The individual metal abundances are listed in TableQland provide [(Xc + Xn + Xq -H XNe)/(XMg + Xsi + Xs + 
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Xca + -'^Fe)] = 0.14 in case of the variable enhancement mixture, and 0.007 for the constant enhancement one. 
The total metallicity ranges from Z = 0.011 to 0.048 . For t he calculation of the latter tables the new low- 
temperature molecular opacity code bv lFerguson et alJ ( l2005l FA05) and for the high-temperature regime the 
OPAL online calculation of opacity tables' were used. We therefore are faced with the fact that the new tables 
(for log T < 4.5; but used only for log T < 4.1) differ from the old ones both in the internal metal distribution 
("mixture effect") and in the method of calculation ("generation effect"). For the higher temperatures only the 
mixture effect is present. 

The purpose of the present paper is to present new low-mass stellar models up to the red giant tip and the 
influence of variations in a-enhancement on the evolution at given total super-solar metallicity. Second, we 
investigate the influence of the new low-temperature opacities and that of a variation of internal element dis- 
tributions in these tables. This is an important issue because of the fact that one cannot expect to have always 
opacities at hand which are fully consistent with the mixture chosen for the stellar models. For completeness, 
we also discuss metal-poor cases. 

The stellar models and the sets of opacity tables are introduced in Sect. 2. The basic results of the cal- 
culations and various tests are the subject of Sect. 3. The low-metallicity case follows in Sect. 5, before the 
discussion and conclusions close the paper in the final section. 

2. Model details 

2.1. Stellar evolution program 

The stellar mod e ls we re calculated with the Ga r ching stellar evolution code, which was described in 
IWeiss & Schlattf flOQ&i and lately in IWeiss et aP ll2005l) . We will concentrate here on the description of 
the opacity tables. Our code uses tables of Rosseland mean opacities k for mixtures quantified by the mass 
fractions of hydrogen and metals (both ranging from to 1 ; in total of order 80), and with a temperature and 
density grid of about 85 and 25 grid points^. The interpolation in this grid is done by a two-dimensional, 
bi-rational spline algorithm (Spaeth 1973). In mixture we use parabolic polynomials first in X (hydrogen) 
between the three tables closest to the actual value, and then in logZ (metallicity). In practice, and as long as 
the total metallicity is not changing, e.g. either by advanced nuclear burning phases or metal diffusion, tables 
for only three metallicity values are sufficient for the whole calculation. For evolutionary stages from core 
helium burning on, we have special core tables, but these are not of interest for the present work. 

The tables themselves are the end product of the merging of source data. We have four main regions in 
the T - R domain: 

- log T < 3.8 , for which we use the Wichita State Alexander 8c Ferguson molecular opacity tables men- 
tioned in the introduction 

- 8.7 > log r > 4.1, for which OPAL tables a re used 

- high density: here we employ the results bv lltoh et al.l (Il983l) for electron conduction opacities 

- log r > 8.7, for which no OPAL data are available; here we use the old Los Alamos Opacity Library (see 
Sect. 1; in the present case this is irrelevant). 

In between these regions are transitions. Between 3.8 < log T < 4.1, where both Wichita State and OPAL 
data are available, we have a linear transition in log k along with log T from one source to the other. The 
agreement between both tables is excellent and the transition is very smooth. At the high-density edge we add 
radiative (k^) and conductive (k^) opacities according to 

1/k = 1/k, + I/Zfe. 

Since with increasing density Kc is dominating, the radiative contribution can be omitted once > k^. 
However, in particular at log T < 5 the radiative tables end before this situation is reached. In case the gap 
between the end of the radiative opacity table and the density from which on Kc is already lower than the last 
radiative value available is only 1-2 dex, we boldly interpolate over this gap (cubic spline). If the gap is too 
large, the final A--table has to end here. Should we run out of the table definition range during the stellar model 
calculations, we use the last table value. This happens at isolated points in some low-mass main-sequence 
envelopes and cannot be avoided as long as the input tables do not cover the full p - T-plane. 

All our tables are constructed in this same way in a separate step before they are used by the stellar 
evolution code. 

' http : //www-phys . llnl . gov/Research/OPAL^ 

^ As the density coordinate we actually use the usual logR = logp - 3 log T + 18 
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Table 1. Chemical abundances of the metals. Columns 2 and 3 contain the standard fOrevesse & Sauvail ( ll998l) 
solar abundances on a logarithmic scale where the hydrogen abundance is 12.00 and as rela tive mass fractions. 
Columns 3 and 4 show the abundances of the a-enhanced mixture by ISalaris & WeissI (il99&. "a-v") and 
columns 5 and 6 for a constant a-enhancement of 0.4 dex ("a-c"). a-elements are printed in italics. 
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2.2. Stellar models 

For the larger project about stellar population synthesis (Coelho et al. 2006, in preparation) we calculated 
mixtures of varying metallicity from Z = 0.005 to 0.048, but in the present paper we will be concentrating on 
one mixture, which is X = 0.679, Z = 0.032. Al l g-elements ar e enhanced over the solar-scaled abundances by 
0.4 dex, for which we used the solar mixture bv lGrevesse & SauvaU1998i.) . For comparison with a solar-scaled 
mixture with identical [Fe/H] models with the mixture X - 0.735, Z = 0.017 was chosen; the higher total 
metallicity is thus assumed to be connected with an increase in helium abundance. In mass our models cover 
the range from 0.6 to 1.2 Mq in steps of 0.05 M©. Our calculations start with the zero-age main-sequence 
(ZAMS) defined by homogeneous stellar composition and vanishing gravo-thermal energies in the initial 
model. Since this composition is not identical to nuclear equilibrium in the core (^He, CNO-isotopes) the 
initial few million years are spent to adjust the composition accordingly. This leads to a well-known short 
loop around the ZAMS position, which will be visible in some of the diagrams. We follow the evolution into 
the core helium flash at the tip of the RGB, until the helium luminosity is of order 1000 Lq. 

We emphasize that the chemical composition in the models is always the same (a-c), as given in Tabled 
What is varied in the calculations presented below, is just the internal composition of the opacity tables. 

2.3. Opacity tables 

The fi rst s et of models (A) was calculated with the same opacity tables as those used bv ISalaris & WeissI 
(Il998l) and'Salasn ich et alJ tOOd) . which are a combination of OPAL and Wichita State opacities for the a-v 
mixture^ They were c alculated at the time when the new Wichita State low-temperature molecular opacities 
jFerguson et all2005l) were not yet available, in particular not for non-solar metal compositions, and under the 
assumption that the slightly different table composition would not influence the model properties significantly, 
such that this internal inconsistency could be tolerated. 

Set B followed after low-T tables for mixture a-c had been calculated. High-T OPAL tables were obtained 
from the OPAL-website for the identical mixture. This is the most recent set of opacity tables and fully 

' These tables were originally provided by F. Rogers and D. Alexander by private communication 
(1994/1995), but later made public on the respective websites, http : //www-phys ■ llnl . gov/Research/OPAL and 
|http: //webs. wichita.e du/phy sics/opacity7] 
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Table 2. Details about the tables of Rosseland mean opacities for the various sets of stellar models. Low- 
and high-temperature ranges are as described in Sect. 12. II Mixture labels are as in TableQland "code" entries 
refer to the relevant publication, where the calculation method and details are described. "AF05" refers to a 
re-computation in 2005 using the original AP94 program. 
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consistent with the model compositions. The models for the population synthesis project (Coelho et al. 2006, 
in preparation) have been computed with it. 

To analyse the results of the following section, additional sets of tables were produced for various test 
cases. They consist of inconsistent combinations of low- and high-T tables and are listed in TableElas well. 
The final case R is a repetition of case A, but the tables were recomputed with the same code bv lFerguson et alJ 
f2005.) . 



3. Results 

3.1. From old to new (Case A to B) 

The scientifically interesting question is how much the use of opacity tables with an internal a-element distri- 
bution (a-v) differing from that of the models themselves (a-c) influences the evolution. This is an interesting 
aspect of stellar modelling as long as hghly accurate opacities cannot be calculated on-line with the models. 
To evaluate this, we computed cases A and B, the latter being the fully consistent one. We emphasize again 
that the model composition in case A is the same as in B, but that the internal metal distribution of the opacity 
tables is slightly diff'erent. The resulting tracks are shown in Fig.^ The eff'ect is drastic: for the 0.9 Mq model, 
ZAMS and TO-positions are 0.08 resp. 0.06 dex brighter, the turn-off (TO) also being 0.02 dex cooler and the 
main-sequence lifetime is 3.3 Gyr shorter than that for case B (17.8 Gyr). RGB temperatures are higher by 
* 250 K, too. Theses changes are also typical for all other masses. 

In the tracks, small kinks are visible both on the main sequence and on the subgiant bran ch. The reason 
is that in the phase lying between these points the OPAL equation of state jRogers et al.l 19961) does not cover 
the conditions in the outer stellar envelope and we switch back to a simple Saha-type equation of state. This 
effect is less pronounced for lower stellar masses. We have repeated a few cases using the new OPAL-based 
equation of state by Irwin (described in ICassisi et alJl2003l) . which is more extended. The tracks from this 
repeated calculation agree extremely well with those shown, but in addition are smooth without any kinks. 

Since the strong influence of the different ff-enriched opacity tables could be due to either the generation 
or the composition effect, or both, we set out to isolate this by additional tests in which we used specific 
opacity tables with only one effect being present. 

3.2. T1 : new low-T opacities 

For this test case we used only the low-T opacity tables for the a-c mixture. The high-T tables are still those for 
mixture a-v as in case A. Figure |5] shows the resulting tracks for three different masses. Temperatures along 
the RGB (for M = 0.8 and 1.0 Mq) are again lower by more than 200 K for the new low-T tables, while the 
ZAMS and TO luminosities differ only by a few hundreths of a magnitude. Correspondingly, main-sequence 
lifetimes agree within 0.4 Gyr (0.8 Mq and 1.0 Mq). The new low-T opacities apparently affect the envelope 
structure significantly, while the high-T tables with the new mixture do influence the core properties. 
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Fig. 1. Comparison of evolutionary tracks for cases A and B. Masses range from 0.6 Mq to 1.2 Mq in steps of 
0.1 Mq (we omitted the intermediate mass values). The calculations were stopped either at the RGB tip or if 
the age exceeded 50 Gyr For this reason the tracks do not extend to the RGB for M = 0.60 Mq for case A, 
and M - 0.60 and 0.70 for case B. The composition of both model sets is identical (a-c), but that of 
the opacity tables is not. 




3.8 3.7 3.6 3.5 3.4 

log Teff 



Fig. 2. Comparison of evolutionary tracks for cases A and Tl . Masses are 0.6, 0.8 and 1 .0 M©. 



3.3. T2: new high-T mixture 

Since the method for calculating the OPAL opacities has not changed, replacing the high-T tables with those 
for the new mixture does not imply a generation effect. Again, we show the comparison between case A and 
this test case (T2) in Fig.|3 Since the low-T opacities are now identical, the RGB temperatures completely 
agree. However, changes on the main-sequence remain. The ZAMS luminosity is lower by 0.07 dex in case 
T2 (M - 0.80 Mq), and the MS lifetime longer by 4.3 Gyr. The reason for this effect is evident from Fig.|3 
where we show a contour plot of the ratio of Rosseland opacities between the old and new composition. Over 
a wide temperature and density range, typical for core temperatures of main sequence stars the a-c opacities 
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Fig. 3. Comparison of evolutionary tracks for cases A and T2. Masses are 0.6, 0.8 and 1 .0 M( 




Fig. 4. Ratio of Rosseland opacities k between a typical table used in cases A and T2. The composition is 
X - 0.700, Z = 0.03. Shown are isocontours of the ratio. 



are higher by 5-20%. The differences are very similar also for lower helium content. An increase in opacity, 
as is well known, leads to a lower luminosity and thus a longer burning time. 



3.4. T3: new generation of low-T tables 

After we identified the main influence of low- and high-T tables, when going from case A to B, we are sep- 
arating generation- from mixture efi ^ect for the low- t emper ature opacities. To this end, tables for mixture a-v 
were calculated using the method of lFerguson et alJ (l2005l) and combined with the a-c OPAL tables. We thus 
compare to case B, with only the mixture-effect at low temperatures present, shown in Fig. [S] Obviously, 
in spite of employing mixture a-v for the molecular opacities, the track resembles that of case B, the com- 
pletely consistent a-c calculation, very much. In fact, for M = 0.80 and 1 .0 Mq ZAMS luminosities agree 
within 0.003 dex, TO luminosities and effective temperatures better than 0.004 resp. 0.003 dex, and RGB tip 
temperatures better than 0.09 dex. The TO and RGB tip ages are within 100 Myr. The conclusion from this 



8 



Weiss et al.: a-enhanced opacity tables 




experiment and the one of Sect. l3.2l is therefore that for the low-temperature opacities the generation effect is 
the dominating one and that within the same generation of tables the mixture effect is almost negligible. This 
is different from the effect on the atomic opacities of the interior (Sect. l33t . 
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Fig. 6. Comparison of Rosseland mean opacities in a table with mixture X - 0.7, Z - 0.02, taken at log R - -3 
as function of temperature (upper panel). Lower panel: differences in log k. F05 denotes the new calculation 
of tables for the same a - v element ratios as in the original AF94 table. 



3.5. R: low-T tables recalculated 

The result of Sect. 13.41 is quite surprising given the fact that iFerguson et alJ (l2005l) reported rather small 
changes in the new molecular opacities for solar-type mixtures. To verify this conclusion further, we in- 
vestigated the low-temperature tables for a-\ mixtures in more detail. A comparison displayed quite large 
differences of ±10% around \ogT x 3.6, increasing further to lower temperatures. Fig. |6l shows a sample 
comparison at log R = -3 for the mixture X = 0.7, Z = 0.02 between the old AF94 and the new FA05 table. 
For log r < 3.4 the difference in fact increases up to a factor of 3. This finding emphasizes that the new cal- 
culation programs alone are responsible for the change in opacities. Using these recalculated FA05 tables for 
mixture a-v, combined with the corresponding high-T OPAL tables, we show the comparison of this case R, 
which is the "2005 update" of A, with case B in Fig.Q We thus see here the effect of changing the a-element 
enhancement factors within the opacity tables, which, however, were calculated with the same program in 
both cases. Obviously, along the RGB the two cases agree very closely, while the differences in the earlier 
evolutionary phases is that of Sect. l3.3l and Fig.|3l 

3.6. a-enhanced low-T tables of 1994 

In contrast to the solar metallicity case investigated in iFerguson et al.l ( l2005l) . the re-calculation of the o'-v 
low-temperature tables with the new molecular opacity code revealed the large changes of Fig.|6j in particular 
at r < 2500 K. The original calculations by D. Alexander (1994, pr i vate co mmunication to A. Weiss) were 
therefore repeated with the original code of [Alexander & FergusonI d 19941) . The result is also displayed in 
Fig-Elas the dashed line, and differs strongly from the original 1994 result, shown as the dotted line, in the 
temperature region between 3.25 < log T < 3.35, while it agrees very well with the result from the new code, 
at least down to log T 3.15, below which the differences can be explained by by updated grain physics, e.g. 
more grain species in the equation of state and updated grain opacities, in FA05. 

The conclusion is that an error was made in the original calculations provided to A. Weiss in 1994. This 
error was essentially a series of typos in the input abundances to the opacity tables run at that ti me. The a- 
enhanced tables used in several publications (Salads & Weiss 1998; Weiss & Schlattl 2000: Salasnich et all 
\200(i and others) are therefore erroneous and should not be used any longer. At least at solar and super-solar 
metallicities they influence the effective temperatures of the models strongly by the underestimation of the 
opacity at low temperatures: Red Giant models are rendered too blue . Similar opacity tables fo r dif ferent a- 
enhancement factors calculated with the same code and used, e.g. in IVandenBerg et all ( l2000l) and I Yi et"al] 
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do not suffer from this error We also emphasize that most of the papers using the erroneous tables 
dealt with significantly lower total metallicity at which the models are much less sensitive (see below). 
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4. Low-metallicity models 

Figure|Slshows the same opacity comparison between a-c ("f05") and a-v tables, the latter for the new ("af94 
corr") and original ("af94 old") calculation, for Z = 0.0001. The situation resembles that of Fig.|6l but the 
opacity itself is lower by an order of magnitude. All changes should therefore result in milder effects on the 
stellar models. 




Fig. 8. As Fig.il butforZ = 0.0001. 



Indeed, Fig. |9] (lower panel) demonstrates this for the case of M = 0.8 Mq, and Z - 0.001: The track 
labelled " AF94" wa s calculated with the original a-enhanced table and corresponds to the tracks used in, 
e.g. Sa laris & Weissl Q298). The same code as in this paper was used, too. Employing the newly calculated 
opacities (F05) again results in a cooler RGB; all other properties remain virtually unchanged. However, the 
difference in Teff amounts to less than 150 K compared to about 250 K for the Z = 0.032 cases discussed so 
far 

The cooler RGB temperatures imply t hat the fits to real globular cl uster data will have to be reconsidered 
(they were, in fact, fitting quite well; see'Salaris & Weiss" 1997', '1998'), a deeper look into this is war ranted. 
With the new generation of low-T tables also the opacities for solar mixtures changed slightly ( Ferguson et al.l 
l2005h . which implies that the solar models need to be recalibrated and the parameter cmlt in the mixing 
length theory needs some adjustment. Using the Garching stellar evolution code A. Serenelli (private com- 
munication, 2006) found that omlt had to be increased from 1.78 to 1.88. Similarly, for the code used for 
the models of Fig.|9l a solar calibration (without diffusion) resulted in an increase by 0.12 to the value of 
ffMLT - 1-95. Using this value, the dotted track in Fig. ^results. This would be the completely consistent, 
correctly calculated track to be used in an update of globular cluster age determinations. Obviously, the new 
opacity tables along with the consistent aMLT is very close to the AF94-track, and the effect on isochrones 
is minimal (e.g. effective temperatures differ by 30 K at logL/L^Mw - 2)^ whi ch is a very convenient, but 
fortunate fact. Thus the globular ages determine in ISalaris & Weissl ( Il998l) and ISalaris & We iss (2002) are 
correct and remain unchanged. We also investigated how a recalibrated q;mlt would modify the effect on the 
more metal-rich effect. This will be reported below and in Fig. 

The second aspect to be considered is that of the mixture effect in the high-T opacities (Sect. l33t . which 
could affect age determinations of globular clusters by the increased main-sequence lifetimes. The difference 
in the opacities, going from one a-enhanced mixture to the other, is vanishing with decreasing metallicity 
(Fig.llO>. The difference for the lowest metallicity considered here amounts to less than 3% everywhere. As a 
consequence, the turn-off ages agree within 200 Myr, or about 2% of the main-sequence lifetime. 

The evolution shown in the upper panel of Fig.|9ldemonstrates the negligible effect of varying the internal 
element distribution on the HRD of the 0.8 Mq star: at fixed Z = 0.001 the three tracks using opacity tables 
for solar (dash-dotted line), a-\ (long-dashed) and a-c (dotted) element ratios are extremely close to each 
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Fijg. 9. Lower panel: Evolution ary tracks for an M - 0.8 Mp star with Z = 0.001 calculated with the original 
[Alexander & FergusonI ( Il994i AF04) and the new lFerguson et all ( l2005l F05) code for the a-w mixture. The 
dotted line refers to a track with the new opacities, but a recalibrated mixing parameter aMcr- Upper panel: 
Effect of using opacity tables with different internal element distributions but identical Z; see text. 



other, the largest - but still very small - effect appearing at the turn-off: The effective temperature agrees for 
the scaled-solar and a-c mixture within 5 K, while the g- v track is ^ 30K ho tter Luminosities are higher 
by 0.005 resp. 0.014 dex. This confirms again the results bv lSalaris et alJ ( Il993l) . In contrast, the "generation 
effect" of the AF94 track is much larger. No recalibration of q-mlt was done here; this would shift all new 
tracks in the same way as the dotted one in the lower panel. 

Given these results, we return briefly to the super-solar case. In Fig. ^2 we compare various cases for 
two mass values (0.8 and 1.2 Mq). The solid line corresponds again to case B (a-c mixture both in the model 
and in the opacity tabl es). The fully equiv alent models for the a-v mixture, taken again into account in both 
model and opacities ( F erguson et alJ20 05') is indicated by the dotted line. In both these cases ckmlt = 1.6. The 
recalibration of the mixing length parameter yields a value close t o 1.7; this ca se is represented by the dashed 
line. Finally, case A (model composition a-c, old Alexander & Fer gusonI 1 994l tables, otmlt - 1 -6) is the most 
extreme case (dash-dotted line). In contrast to the low-metallicity case, the recalibration of q-mlt does not 
approximate case A. It is also obvious that the change of mixture has its largest effect on the main-sequence 
due to the increased high-temperature opacities, but less so on along the RGB, where temperatures are much 
more affected by q-mlt, which in turn has a smaller effect along the main-sequence, in particular of course for 
the higher mass with very thin exterior convection zone. 





Fig. 11. Overview over various cases for two mass values (0.8 and 1.2 Mq) for our standard composition of 
X - 0.679, Z - 0.032. See text for more explanations. 
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5. Discussion 

The diversity of chemical compositions for stellar models will remain to be much larger than that of available 
opacity tables, which, however are crucial for accurate model calculations. We have therefore investigated 
how opacities calculated for similar, but not identical a-element enhanced mixtures influence the evolution 
of stellar models of low mass and super-solar metallicity, where effects are expected to be large. While the 
models had constant a-enhancement factors (a-c; Tabled we used tables for the same composition, but also 
for another one of varying enhancement factors, a-v. The latter opacity tables had been computed about a 
decade ago with the [Alexan der & Ferguson ( 1994) code in response to a specific request. 

The initial comparison (Fig.[0 revealed drastic changes in the temperature of red giants and the location 
and duration of the main sequence. These changes were much larger than what could be expected on the basis 
of the experience with using solar-scaled mixture opacities instead of or-enhanced ones cSala ris et al. 1993; 
LSalaris & Weissll998l) . 

The origin of low- and high-temperature opacities dififering, we selectively replaced them to identify the 
cause for the diff'erences found. The lower stellar temperatures found in case of a- c opacity tables are, i n 
particular along the RGB, almost exclusively due to the low-T molecular opacities Jperguson et al.lE005l) . 
while the main sequence luminosity and thus lifetime varies with the composition of the (OPAL-)high-T 
tables, the latter by up to 20% ! As the code computing the latter tables has not changed since Iglesias & Ro ger j 
(Il996l) . this was identified as a pure composition effect: The a-c opacities, for given (X, Z), is consistently 
higher in a temperature and density range covering the whole energy producing core of low-mass main- 
sequence stars, which explains the lower luminosity found. We have performed various tests to make sure that 
we obtained and used the correct tables and have thus verified that the individual a-element abundances play 
a significant role for the opacities at stellar core temperatures. In addition, we created opacity tables for the 
same mixtures using the latest Opacity Project data (Badnell et al. 2005; Seaton 2005). While in the interesting 
temperature range the OP Rosseland mean opacities differ from the OPAL ones, the difi^erential effect going 
from a-v to a-c abundance ratios is very similar. If the opacities we have been using are confirmed to be true, 
the conclusion is far-reaching: at solar- and super-solar metallicities, the lifetime of ff-enriched low-mass 
stars is sensitive to the individual abundance pattern. Accurate age determinations of stars or stellar systems, 
such as elliptical galaxies or of the galactic bulge, may require very detailed knowledge about the chemical 
composition. 

For the low-T opacities we further had to discri minate between the com position and the generation effect, 
as the a-c tables had been computed with the recent 'Ferguso n et alJ (12005^ code. As the result of our tests we 
discovered that the original a-v tables suffer from an error and are incorrect. Recomputations using both the 
new and the old Wichita State molecular opacity code resulted in very similar opacities, which, when used 
(Sects. I3!4l and l3.5i . yield much lower RGB temperatures."* 

The conclusion of our investigation is that at high metallicities both the error in the old low-T tables and the 
composition details in the high-T opacities significantly influence our models. Previous calculations using the 
erroneous tables at solar metallicities should be repeated with the new tables for the same a-v composition. 
For sub-solar metallicities both effects tend to decrease, and in particular for typical Pop. II compositions 
variations in individual element abundances are no longer significant for lifetimes. Age determinations of 
globular clusters are thus ind ependent of composition details. We therefore confirm once again the basic 
results of ISalaris et alJ ( Il993h . The erroneous low-T tables, in contrast, still render the tracks too hot on the 
RGB, such that colour-based age determinations would be affected (yielding ages too low). Incidentally, 
however, a recalibration of t he mixing le ngth parameter aMu, which may be considered necessary when 
using the new low-T tables of lFerguson et aL (.2005) almost recovers the old RGB temperatures (Fig. 13 lower 
panel). 

To summarize, we found that metal-rich low-mass stellar models are very sensitive to details of the chem- 
ical composition used in opacity tables. It is therefore advisable to employ opacity tables consistent with the 
assumed stellar chemical composition. The range of massive stars still needs to be investigated. For popula- 
tion synthesis models of metal-rich galaxies knowledge about the mean ff-enrichment may not be sufficient 
for any reasonable statements about the star formation history. 

Acknowledgements. JF acknowledges support from NSF grant AST-0239590 with matching support from the State of 
Kansas and the Wichita State University High Performance Computing Center funded by NSF grants EIA-0216178 EPS- 
0236913. 



The old opacity tables for this mixture have been removed from the web-site 
|http : //webs . w ichita . edu/phy sics/opacity/| and replaced by new calculations. 
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